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Q , Abstract A recurrence equation is a discrete integrable equation whose solu- 

tions are all periodic and the period is fixed. We show that infinitely many 
CN) ' recurrence equations can be derived from the information about invariant va- 

rieties of periodic points of higher dimensional integrable maps. 
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J© ' 1 Introduction 
o 

<N ■ A recurrence equation is a discrete integrable equation whose solutions are all periodic 
\q \ and the period is fixed. Some of them had been known for some years, while some others 
have been found recently. In this contribution we would like to show that infinitely many 
recurrence equations can be derived from the information about invariant varieties of 
periodic points of higher dimensional integrable maps. Especially recurrence equations 
associated with the Quispel, Roberts and Thompson (QRT) map[lj are shown to exist 
one for each period. 

Some examples of the recurrence equations are [2] 



Oh 

-i— > 



^ ; x n+ i = — , a : constant (1) 

l + X 

x n+l — > \ z ) 
Xn—1 

1 -|- X n -\- X n —\ , . 

X n +1 = • (3) 

An interesting feature of these equations is that, for an arbitrary initial value, the solution 
is always periodic with a fixed period. The period is 2 in the case of ((H), 5 in the case of 
and 8 in the case of These equations were named the recurrence equations by the 
authors of [31 H] who found many other examples of this type recently. 

Apparently the recurrence equations are integrable. However there has not been 
known, to our knowledge, any method to find them systematically. The purpose of this 
article is to develop a method to derive the recurrence equations from integrable maps in 
a systematic way. 

Our key observation is that, writing (x n , x n _i, x„_ 2 ) as (x,y,z), the recurrence equa- 
tions (pQ), (|2J), ([3]) are equivalent to the higher dimensional maps 



(4) 



(x,y) -> (X,Y)=^±^, (5) 
(x,y,z) -> (X,y,Z)= + £ + x> y V ( 6 ) 



respectively. There are a pair of fixed points at x — ±\/a, a; = y — (1 ± \/5)/2, and 
x = y = z = 1 ± \/2 for each map (BJ, (jSJ), (JS). Otherwise an arbitrary point on the 
complex space can be an initial point of the periodic map of the corresponding period. 

We have shown, in our recent paper (5j [6], that periodic points of higher dimensional 
integrable maps with some invariants form an invariant variety for each period. The 
invariant variety of periodic points is determined by imposing certain relations among the 
invariants. Every point on an invariant variety can be an initial point of the periodic map 
of the same period. All images of the map stay on this invariant variety. Therefore the 
map defines a recurrence equation of the fixed period if it is constrained on the invariant 
variety. In some cases the invariant varieties can be derived iteratively for all periods. 
We can associate one recurrence equation to every invariant variety, thus obtain infinitely 
many recurrence equations. 

We explain briefly the notion of invariant varieties of periodic points in §2. Many 
recurrence equations associated with the invariant varieties will be derived in §3. We 
discuss, in §4, a method which enables us to derive series of recurrence equations. 

2 Invariant varieties of periodic points 

Let us consider an iteration of a rational map on C d , where C = {C, oo}, 

x = (x 1 ,x 2 ,...,x d ) -> X=(X 1 ,X 2 ,...,X d )=:X«, (7) 

and assume ifi(x), /^(x), H p (x) be the p invariants. We are interested in the behaviour 
of periodic points satisfying the conditions 

X (n) =x, n = 2,3,.... (8) 

If hi, h 2 , • • • , hp are the values of the invariants determined by the initial point, the 
map ([7]) is constrained on the d — p dimensional algebraic variety V(h), 

V{h) = {x| Hi{x) = hi, i = 1,2, ...,?}, (9) 

and the periodicity conditions (jHJ) are reduced to the constraints on some d — p functions 

p(n). 

^{h, h 2 , hp, £i, 6, = 0, a = 1, 2, .., d -p, n>2. (10) 

Here by £1,^2, ■■■,£d-p we denote the variables which parameterize the variety V(h) after 
the elimination of the p components of x. Note that the fixed point conditions (n = 1) 
are excluded in fflOl) since they are nothing to do with the invariants. 

For an arbitrary set of values of hi, h 2 , h p , the functions T^'(h,^) define an affine 
variety, which we denote by V^")((r)), i.e., 

^«r» = {e|rW(/i ) = 0, a = l,2,...,d-p}, n>2. 



In general this variety consists of a finite number of isolated points on V(h), hence zero 
dimension, corresponding to the solutions to the d — p algebraic equations ({TO]) for the 
d — p variables £1,^2, ■■■,^d-p- In this case we say that the periodicity conditions (jHJ) are 
'uncorrelated'. If the values of the invariants are changed continuously these points move 
all together and form a subvariety of dimension p in C d . Needless to say this case includes 
a map with no invariant. 

There are possibilities that the equations ffTUj) impose relations on hi, h 2 , h p instead 
of fixing all £ a 's. Let I be the number of such equations. We write them as 

^(h 1 ,h 2 ,...,h p ) = 0, a = 1,2,..., I, (11) 

instead of to emphasize independence from ^'s. If m is the number of the rest of 
the equations 

T^(h u h 2 ,...,h p ^ 1 ,^,...^ d _ p ) = 0, a = l,2,...,m (12) 

d — p — in variables are not determined from the periodicity conditions. This means that 
V"( n )((r)) forms a subvariety of dimension d — p — m of V(h). We say that the periodicity 
conditions are 'correlated' in this case. In [5] we have proved the following lemma: 

Lemma[5j 

A set of correlated periodicity conditions satisfying min{p, d — p} > I + m and a set of 
uncorrelated periodicity conditions of a different period do not exist in one map simulta- 
neously. 

When m = 0, in particular, the periodicity conditions determine none of the vari- 
ables £i,£2j ■■■i^d-p but impose I relations among the invariants. Then the affine variety 
VA n )((r)) coincides with V(h). In other words every point on V(h) is a periodic point of 
period n, while V(h) itself is constrained by the relations among the invariants. We say 
the periodicity conditions are 'fully correlated' in this particular case. If we replace hi 
by Hi(x) in ^^(h) the periodicity conditions (iTTj) enable us to consider the constraints 
on the invariants as constraints on the variables x. We denote by v^^)) the affine 
variety generated by the functions 7 ( ^ n )(if 1 (x), H 2 (x), H p (x)), and distinguish it from 
V {n) ((T}). Namely we define 

7 W(if 1 (x),^ 2 (x),...,i/ p (x)) = 0, a = 1,2,...,/}. (13) 

We call f^ n ^((7)) 'an invariant variety of periodic points', whose properties can be 
summarized as follows: 

• The dimension of v ^((7)) is d — I (> p). 

• Every point on f ( - n ^(( 7 )) can be an initial point of the periodic map of period n. 

• All images of the periodic map starting from a point of f^™- ) (( 7 )) stay on it. 

• i> ( ' n )(( 7 )) is determined by the invariants of the map alone. 



« (n) «7» = {x 



If the periodicity conditions of period n are fully correlated, i.e., when m = in (fl2l) . 
the condition min{p, d — £>} > Z + m is always satisfied as long as (|lip has solutions. Our 
theorem thus follows from the lemma and this fact immediately. 

Theorem [3] 

// there is an invariant variety of periodic points of some period, there is no set of 
isolated periodic points of other period in the map. 

This theorem tells us nothing about the integrability of a map. To proceed further we 
assume that a nonintegrable map has at least one set of uncorrelated periodicity condi- 
tions. This is certainly true if the map has a Julia set. Once we adopt this observation 
as a working hypothesis, our theorem is equivalent to the following statement: 

If a map has an invariant variety of periodic points of some period, it is integrable. 

In order to support this proposition we have investigated various known integrable 
maps and found invariant varieties of periodic points in all cases if there are invariants. 

3 Derivation of recurrence equations 

If there is an invariant variety of periodic points of period n, every point on the variety 
can be an initial point of an n periodic map. All images of the map are on the variety 
before the map returns to the initial point. Therefore this variety is clearly distinguished 
from the rest of C d and is reserved only for the maps of period n. In other words if the 
initial point is on this variety the map is always n period. 

This fact enables us to derive a recurrence equation once an invariant variety of periodic 
points is known. Let 

x i ~> Xj = f j (x 1 ,x 2 ,...,x d ), j = l,2,...,d (14) 

be the map and (fl3l) be the invariant variety of period n of this map. We solve 

7 i n) (#i(x), # 2 (x), tf p (x)) = 0, a = 1, 2, I. (15) 

for I variables, say Xd-i+i, • and substitute them into fj, j = 1,2, ...,d — I of (1I41) . 
Then every initial point of the map Xj = fj, j = 1, 2, d — I is constrained on v^({'j)), 
thus we obtain a recurrence equation of period n. 

The simplest method to achieve this program is to find the Zth elimination ideal of the 
functions {Xj — fj, j = 1, 2, d — 1} and {yffi, a = 1, 2, I}. If the ideal is generated 
by the functions F^'s satisfying 

F^ n) {X 1 ,X 2 ,...,X d _ l ,x 1 ,x 2 ,...,x d . l ) = 0, j = 1,2,..., d- 1, (16) 

the recurrence equations are obtained by solving (TIB"]) for X 1; X 2 , ...,X d _i. Generally the 
solutions are not rational. But it does not cause any trouble, since the integrability of 
the map has been guaranteed from the begining. For all initial values xi,x 2 , ...,x d -i, the 
solutions of (TTE1) are periodic, hence are integrable. 



For an illustration let us consider the map 



(x,y)^(X,Y) 



xy, 



yi} + x) 
1 + xy 



(17) 



This map has one invariant H(x,y) = y(l + x) and the invariant variety of period 3 is 
given by the zeros of 



7 (3) (x,y) = H 2 + H + l 

= x 2 y 2 + 2xy 2 + y 2 + xy + y + 1. 



(18) 



The 1st elimination ideal of the functions X — xy and ( fTB"l) is generated by the function 

F (3) (X, x) = (x+ l) 2 X 2 + x(x + 1)X + x 2 , 
from which we obtain two maps: 



x 



10- 



x -> X = < 



The iteration of the first map yields 



x + r 

X 

X + 1 



(c 3 = 1). 



X 



UJ 2 X 



X — ► UJ- 



X + 1 



while the second map yields 



x 



X 

x+1 



—L0 2 X + 1 



LUX 



-UJX + 1 



x. 



X. 



In the rest of this section we would like to present various type of recurrence equations 
associated with invariant varieties of some integrable maps. 



The d dimensional Lotka-Volterra map is obtained by solving [7] 
Xj(l - Xj-x) = Xj(l - x j+ i), j = 1, 2, d 



(19) 



for X = (Xi, X 2 , Xd) under the conditions Xj + d = Xj (j = 1,2, ...,d). The invariants 
of this map are given by 

Hk = ^ji,j3,...Jk X 3i x h ' ' ' X jk(^- ~ a; ii-i)(l — x j2-i) ' • • (1 — x ifc-l) 

(fe = l,2,...,[d/2]) (20) 

r =x 1 x 2 ---x d 

Here the prime in the summation J2' of fl20|) means that the summation must be taken over 
all possible combinations ji,j 2 , ■■■,jk but excluding direct neighbours. The total number 
of the invariants is p = [d/2] + 1, where [d/2] = d/2 if d is even and [d/2] = (d — l)/2 if d 
is odd. 

The invariant varieties have been derived in the cases of d = 3,4 and 5 for some 
periods, explicitly|5J. In all examples the dimension of the invariant varieties is p. Hence 
the dimension of the recurrence equations is also p. 



The 3 dimensional Lotka-Volterra map is given by, writing (xi,X2,xs) = (x,y,z) 



1-y + yz 1 - z + zx 1 - x + xy 
X = x- ■ , Y = y- ■ , Z = z- ■ , 21 

1 — z + zx 1 — x + xy 1 — y + yz 

after solving ffl9l) for X, Y, Z. There are two invariants 

r = xyz, s = (1 -x)(l -y)(l - z). (22) 

The invariant varieties of periodic points have dimension 2 and are generated by the 
functions: 

7 (2) = s + 1 

^(3) _ r 2_j_ s 2_ rs _|_ r _|_ s _|_^ 

^( 4 ) = r 3 s + s 3 — 3rs 2 + 6r 2 s + 3rs — r 3 + s 

7 (5) = r 3 s 4 _ r 3 s 2 _ g r 4 s 5 + ^3^6 + 3^ + ^6 + ^5 + 3^4 _ 3^3 

-6rV - r 6 s 3 + 3rV + s 4 + 21s 4 r 2 + 6s 4 r + r 3 s 7 + s 7 
+27s 5 r 2 - 3s 6 r - rV + 21rV - 10r 3 s 3 - Qrs 7 + s 8 

i (23) 

for the period 2,3,4,5,..., respectively. 

From these data we can derive a set of recurrence equations for each period. For the 
period 2 case we find 

if } = {x- 1)X - x, if } ={y- 1)Y - y, 
and the map is simply given as 

\a;-l 2/ — iy 

In the period 3 case we obtain 

if } = (x 2 - x + 1)X 2 + x(xy -2x + y + l)X + x 2 {y 2 -y + 1) 
if } = ((3x 2 - 3x + l)y 2 - {3x 2 - 5x + 2)y + (x - 1) 2 )y 2 

- y((3x 2 -2x + l)y - (2x - l)(x - 1))Y + y 2 (x 2 -x + 1). 



Since the solutions of (if } = 0, F 2 (3) = 0) are two folds the map has two routes: 

x(y + u 2 ) (1 — uj 2 )y(x + u 2 ) \ 



(x,y) -> [uj : 



x + u 2 ' 3xy + (x + y — l)u 2 J 

( (1 - u 2 )x(y + u 2 ) y(x + u 2 )\ 

n w ■ — 5- -> (a;,j/) 



\3:q/ + (x + y — l)ccJ 2 ' y + lu 2 
2 x(y + uj) (1 - u)y(x + u) 



x + lo ' 3xy + (x + y — l)u J 

(1 - u)x(y + u) 2 y(x + uj)\ 

w ; -> (x,y) 



3xy + (x + y — l)u J y + 00 



Similarly we can derive recurrence equations for larger periods, but their complicated 
expressions are not worth to be presented here for our purpose of this paper. 



The 4d Lotka-Volterra map (x, y, z, u) — > (X, Y, Z, U) has three invariants. The 
invariant variety is generated by the function 

j( 2 > = j{ 1 — 2 = x + y + z + u — xy — yz — zu — ux — 2 

in the period 2 case, from which we derive the recurrence equation: 

f[ 2) = (l-x-z)X + x, 
FP = Y-y(l-x-z), 
F 3 (2) = (1-x- z)Z + z. 

This provides an example of three dimensional map of period 2: 

(x,y,z) ( — — -, y(l-x-z), — — -) (x,y,z). 
\x + z — 1 X + z — 1/ 



The N point Toda map is known equivalent to the d = 2N dimensional Lotka-Volterra 
map [7j. In the case N = 3, the map (x, y, z, u, v, w) — > (X, Y, Z, U, V, W) is defined by 

^ zu + zx + wu TA xv + xy + uv yw + yz + vw 

X = y , Y = z , Z = x , 

yw + yz + vw zu + zx + wu xv + xy + uv 

yw + yz + vw Tr zu + zx + wu TTr xv + xy + uv 
U = u , V = v , W = w- 



zu + zx + wu xv + xy + uv yw + yz + vw 

Since this map has four invariants, 

ti = x + y + z + u + v + w, 

h = xy + yz + zx + uv + vw + wu + xv + yw + zu, 
h = xyz, 
t' 3 = uvw. 

the recurrence equations are expected to be four dimensional. The invariant variety of 
period 3 is given by the intersection of the functions [5] 

7i — C l) 72 — 1 2- 

From this data we derive a four dimensional recurrence equation: 

Fi 3) = (x + y + u)X + (u + v + y)y, 

F^ = (u + v + y)Y — (u + v + y) 2 — (x — v)u, 

FP = ,(u + v + y)U + (x + y + u)u, 

F 4 (3) = (v - x)V + (u + v + y)v, 



or writing the solution explicitly, we find the map 



where 



y 

z 

\u } 



( u + v + y \ 

-y— — — 

x + y + u 
<f>(x,y,u,v) 

x + y + u 

—u 

u + v + y 

u + v + y 

v 

x — V 

(j)(x,y,u,v) 



/ x + y + u \ 

' (f>{x,y,u,v) \ 

x — V 



x- 



u- 



u + v + y 
x — v 

x — v 
u + v + y 

x + y + u 

j 

x — V 



(u + v + y) 2 + (x — v)u 
u + v + y 



( x \ 

y 

z 

V u J 



As the last example we consider the Euler top. Let (x, y, z) be the three components 
of the angular velocity of the Euler top. Then the map (x, y, z) — > (X, Y, Z) satisfying 



X = a(Yz + Zy), Y = f3{Zx + Xz), Z = 1 (Xy + Yx) 



(24) 



defines the discrete analog of the Euler top [8], if the parameters (a, (3, 7) are related to 
the three moments of inertia /, J, K of the top by 



J-K 



a 



K -I 



21 ' ' 2J 
This map has two invariants [U [8], [9] 

Ix 2 + Jy 2 + Kz 2 

-"1 = ; n 1 -"2 - 



7 



I — J 
2K ' 



2^2 



Fx 2 + J 2 y 2 + K 2 z 



1 — I3^x 2 ' 1 — (3^/x 2 

from which we have found an invariant variety of periodic points [6] 



,(3) 



3 + 7 — (3- 



IJ KI V 2JK 

which is generated by the function 

7 ( 3 ) = (1 + /?7x 2 + -fay 2 + a[3z 2 ) 2 - 4af3>y(ay 2 z 2 + [3z 2 x 2 + ^x 2 y 2 ) - 4 
in the period 3 case. The recurrence equation of period 3 is then obtained as follows: 



(3) 



,(3) 



where 



- cm/ 2 ) (cm/ 2 + (3jx 2 - 2 - 2q)X - (1 - a-/y 2 + q)x) 
—ay (1 — osyy + q) (a^/y + p^jx — 1 — 2q) 
a((l - [3jx 2 )(ajy 2 + pjx 2 - 2 - 2q)Y - (1 - (3^x 2 + q)y)' 
-px 2 (l - pjx 2 + q) 2 (a~fy 2 + ^x 2 - 1 - 2q) 

q = J(l - cry?/ 2 )(l - (3^x 2 ). 



The map has two routes, 

corresponding to the zeros of (F^\ F^)\ 

1 - a^y 2 + q x ± y^Jaf3{a^y 2 + /37X 2 - 1 - 2g) 

1 — a7|/ 2 o;7|/ 2 + f3^x 2 — 2 — 2q 

1 - /?7X 2 + q V ± x^/ap(a-fy 2 + (3-fX 2 - 1 - 2q) 

1 — /97X 2 a^y 2 + /37a; 2 — 2 — 2q 

The two routes correspond to the forward and the backward maps starting from the same 
initial point. This means that the discrete Euler top can not start its 3 period motion 
unless the direction of the motion is informed. 



4 Series of recurrence equations 

Let /i(x),/ 2 (x),...,/ d _ p (x), i?i(x),...,if p (x) be some functions of x = (x 1 ,x 2 , ...,x d ). If 
<7j(x), j — d — p + 1, d are the solutions of 

Hi(fi, f 2 , ■-, fd-p, 9d- P +i, 9d) = H i (x 1 ,x 2 ,...,x d ), i = l,2,...,p, 

they define a map 

x -> X= (/i(x),/ 2 (x),...,/ d _p(x),ff d _p+i(x),...,^ <i (x)), 

in which i?! (x) , (x) , ...,if p (x) are invariant. 

For example, if /i(x) is a function of Hi,H 2 , # d -i, and /1 is given by 

/i(x) = - Xi), 

we obtain higher dimensional nonintegrable maps which reduce to the logistic map upon 
the elimination of x 2 ,x%, ...,x d by using the invariants. This idea enables us to consider 
many higher dimensional maps all together, just by studying a simple lower dimensional 
one. 



4.1 Mobius map series 

Following to the above prescription we can derive higher dimensional integrable maps 
which reduce to the Mobius map, if a, b, h are some functions of the invariants H 1 , H 2 , H d _ 
and we define f\ by 

The iteration of this map does not change the form of the map but only changes the func- 
tions a, b, h. Since we have assumed that these functions are dependent on the invariants 
i?i(x) alone, the initial values (a, b, h) remain constant through the iteration. If we write 

X {n) = h (n) X + oW 

1 _l_ h(n)rr v ' 



after n steps, the (n + l)th parameters are related to the nth ones by 



>+i) 



a + a^h^ 



= h 



h (n) + ab (n) 



h(n) + a ft(n) ' i + bh( n )a( n ) ' 1 + bh( n )a( n ) ' 

from which we can determine all parameters iteratively as functions of the initial values 
(a, b, h). 

The periodicity conditions of period n for the map fl25|) are now satisfied if the param- 
eters (a, b, h) satisfy 

(a (n+1) , 6 (n+1) , /i (n+1) ) = (a, 6, /i). (27) 

From the construction it is clear that the periodicity conditions do not fix the values of 
the variable x but impose some constranits on the parameters, hence on the invariants. 
Solving (1271) iteratively we find the invariant varieties of periodic points as follows [5J 

;|l + /i = 0} 
:|1 + h + h 2 + abh = 0} 
z\l + h 2 + 2abh = 0} 

z\l + h + h 2 + h 3 + h 4 + abh{3 + (4 + ab)h + 3h 2 ) = 0} 
z\l-h + h 2 + 3abh = 0} 



(2) 


= {x 


(3) 


= {x 


(4) 


= {x 


(5) 


= {x 


(6) 


= {x 



; (28) 

According to our argument in §1 we should have recurrence equations corresponding 
to the invariant varieties (1281) . one for each period. To obtain recurrence equations we 
must specify the invariants of the map in higher dimensions. Although the dimension 
of the map could be chosen arbitrary, we consider here two dimensions for the sake of 
simplicity. The number of the invariants is one in this case. Let H(x,y) be the invariant. 
A two dimensional map, which reduces to (1251) . will be obtained if we fix f\(x,y) and 
H(x,y) as functions of (x,y). For this purpose we further assume simply that a, b are 
constants and the function f\ and the invariant H are given by 

x + a 



fi(x,y) = H(x,y) 



1 + bx 



H(x, y) = y(l + bx). 



Solving H(fi,g) = H(x,y) for g(x,y) we find a map 
(x, y) (X, Y) = I (x + a)y, y- 



bx 



(29) 



1 + by(x + a) J 

This includes ( 1T71) special case. 

Since we have already the information (128|) of the invariant varieties it is not difficult 
to derive a series of recurrence equations associated with the two dimensional map (1291) . 
one for each period, as follows: 

= (l + bx)X + x + a, 

[1 + bx) 2 X 2 + (1 + ab)(l + bx)(x + a)X + (x + a) 2 , 
[1 + bxfX 2 + 2ab(l + bx)(x + a)X + (x + a) 2 , 
[1 + bx) A X 4 + (1 + 3a6)(l + bx) 3 (x + a)X 3 

[I + Aab + aV)(l + bxf{x + a) 2 X 2 + (1 + 3a6)(l + bx){x + afX + (x + a) 4 , 
[I + bxfX 2 - (1 - 3a6)(l + bx)(x + a)X + (x + a) 2 , 



i?(3) 
i?(6) 



+ 



To convince ourselves let us see some of the maps explicitly. The map of period 2 is 
generated by F^ 2 \ from which we find 



x 



x + a 
1 + bx 



x. 



We notice that the generating functions of period 3, 4 and 6 cases are similar. There are 
a pair of routes for each period. The map in the period 3 case, for example, is given by 



x 



x + a 
1 + bx 

x + a 



x + a\i_ 
/U_ + bx 

x + a/i + 



x. 



where 



1 + bx fi + + bx 

(i± = -(l + ab ± y/(3 + ab)(ab - 1) 



4.2 Biquadratic map 

By studying various higher dimensional integrable maps which reduce to a one dimensional 
map x — > X, we found, in [5], that many of them reduce not to the Mobius map but to 
the 'biquadratic map' defined by the equation: 

aX 2 x 2 + b(X + x)Xx + c(X - xf + dXx + e(X + x) + / = 0. (30) 

Here 

q=(a,b,c,d,e,f)eC 6 (31) 

are functions of d— 1 invariants. The function /i(x) is determined by solving (1301 for X. 

Because of the symmetry of the equation fl30l) under the exchange of X and x, the 
iteration of the map leaves the form of the map and changes only the parameters q, as it 
was shown in [5]. After repeating the iteration n times, the (n + l)th parameters q( n+1 ) 
are determined by and the initial values q. Since the parameters are functions of the 
invariants alone, the periodicity conditions q( n+1 ) = q impose some constraints among 
the invariants different for each period. 

Solving the periodicity conditions iteratively we have found [Sj a series of invariant 
varieties of periodic points, one for each period. If v^ n ' = {x| y( n ' = 0} is the invariant 
variety of period n, the generating functions are given by 

7 (3) (q) = af - be - 3c 2 + cd, 

7 (4) (q) = 2acf - adf + b 2 f + ae 2 - 2c 3 + c 2 d - 2bce, 

7 (5) ( q ) = a 3 / 3 + ( - cfd + 2c fe 2 + fde 2 - 3ebf 2 - e 4 - c 2 f 2 ) a 2 

+ ( - 13c 4 / + 18c 3 /d + de 3 b + 2c/ V + 7dc 2 e 2 - ce 2 d 2 - 2ce 3 b 

+ 2c 2 feb - 7fd 2 c 2 - Uc 3 e 2 + cd 3 f + fb 2 e 2 + fdb 2 - ebd 2 f)a 

- cd 2 b 2 f - b 3 e 3 - Ac 3 deb + cdb 2 e 2 + 13ec% - fb A + 7fb 2 c 2 d 

+ c A d 2 - 5c 5 d + 5c 6 - 2/6 3 ec - e 2 c 2 b 2 + eb 3 df - Ufb 2 c 3 , (32) 



and so on. 



In order to derive recurrence equations, we must specify the higher dimensional maps. 
As we have shown in [5l [6] the symmetric version of the QRT map [TJ, the 3d Lotka- 
Volterra map of (!2li) . the discrete Euler top, a special case of the g-Painleve IV map 
belong to this categoly. For example the 3d LV map ( 12T1) is equivalent to the biquadratic 
map if we choose 

a — r + 1, b — s — 2r — 1, c = r — s, 
d = s 2 + rs + 5r — 2s + 1, e = — r(s + 1), / = 0, 

from which we could derive the invariant varieties (1231) . 



4.3 The recurrence equations derived from the QRT map 

In the rest of this section we want to derive recurrence equations generated from the QRT 
map. Let us consider the two dimensional map 



rr _ (XY gW(y) - pW(y) - x(p'(y)?'(y) - ZWiv)) 



Here 



£'(x) := ax 2 + b'x + c', £"(x) 
rf{x) := b'x 2 + (d' - 2c')x + e', rf{x) 
p\x) := c'x 2 + e'x + f, p"(x) 



&V + (<F - 2c")x + e", 
cV + e"x + /", 



and q' = (a', b', c', d', e', /') and q" = (a", 6", c", d", e", /") are constants. If we write 
(Y, y, x) as (a;( n+1 ), a^"), x^ n ~ 1 '), this is nothing but the symmetric case of the well known 
QRT equation [TJ 



i 


>(")) p "( x W) 


-p' 


{x^)r)"(xW) 


- 


V(a;< n >)£" (i (n) ) 




>("))/(^ n) )) 


p 1 


[xV))Z"(x( n )) 






_ x (n-l)( 


^(i(»))?|"(3;W) 




>W)^"(xW)) 



The map (133]) has an invariant [TJ 



£,'(x)y 2 + r]'(x)y + p'(x) 
£,"(x)y 2 + T]"(x)y + p"(x) 



H ^ V) = - nn^.i I 1^737 I IwTZv ( 34 ) 



hence it can be reduced to one dimensional map x — > X. The calculation of the 2nd 
elimination ideal is rather trivial in this case. If y and Y are eliminated by using the 
invariant H(x,y) = h, the result we obtain is 

£(x)X 2 + rj(x)X + p(x) = (35) 

where 

£(x) := ax 2 + bx + c, r)(x) := bx 2 + (d — 2c)x + e, p(x) := cx 2 + ex + f, 

with 

q = q' + hq. 



If we identify q = (a, b, c, d, e, /) with those of ([3"Ii) . the map (|35|) is exactly the biquadratic 
map (1301). 



In the theory of the QRT map the formula ( 1351) is called an invariant curve [T]. The 
problem of solving the equation (1331) is now converted to finding the coefficients of (|35|) 
iteratively. Our general formula (|32|) gives us the explicit expressions of the invariant 
varieties of the symmetric QRT map ( 1331) . Under these circumstances the problem of 
deriving the recurrence equations is quite simple. Namely we are already given in f[3"2"j) 
the generating functions of invariant varieties in the form 

7 W(q' + F(x,y)q") =0, n = 3,4,5,.... (36) 

The elimination of y from the map (1331) and (|36|) amounts to replace ?/ by X in (|36|) . since 
X = y. Thus we have found that the recurrence equations derived from the QRT map 
are generated by the functions 

j^{q' + H{x,X)q") = 0, n = 3,4,5,.... (37) 

The recurrence equation of period 3 is, for example, 

= (a' + H(x, X)a")(f + H(x, X)f") - (b' + H(x, X)b")(e' + H(x, X)e") 
-3(c' + H(x, X)cf + (c + H(x, X)c") {d! + H(x, X)d"). 

For some particular choices of the QRT parameters (q', q"), recurrence equations had 
been derived by the authors of [3 J . Our formula (1371) provides such equations for all QRT 
parameters and for all periods. 

Acknowledgements 

The authors would like to thank the organizers of the SIDE VII meeting held in 
Melbourne, who extended us a kind hospitality during the meeting and gave an oppotunity 
to write this contribution. 

References 



[i] 

[2] 

[3] 
[4] 

[5] 
[6] 



Quispel G R W, Roberts JAG and Thompson C J 1988 Phys. Lett. A 126 419, 
1989 Physca D 34 183 

Graham R L, Knuth D E and Patashnik O 1994 Concrete Mathematics (Addison- 
Wesley) 

Hirota R and Yahagi H 2002 J. Phys. Soc. Jpn. 71 2867 

Hirota R and Takahashi D 2003 Discrete and Ultradiscrete Systems (Tokyo: Kyoritsu 
Shuppan) (in Japanese) 



Saito S and Saitoh N 2006 ArXiv: math-ph/06 10069, to be published in J. Phys. 
Soc. Jpn., 2007, V.76 No.2 



Saito S and Saitoh N 2006 ArXiv: |math-ph/0610083, to be published in Symmetry, 



Integrability and Geometry: Methods and Applications (SIGMA) 



[7] Hirota R, Tsujimoto S and Imai T 1993 in Future Directions of Nonlinear Dynamics 
in Physical and Biological Systems, ed. by P.L.Christiansen et al., (New York: Plenum 
Press) p 7 

Hirota R and Tsujimoto S 1995 J. Phys. Soc. Jpn. 64 3125 
[8] Bobenko A I, Lorbeer B and Suris Yu B 1988 J. Math. Phys. A: Math. Gen. 6668 
[9] Hirota R and Kimura K 2000 J. Phys. Soc. Jpn. 627 



